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Abstract 

Polymer  electrolyte  fuel  cell  (PEFC)  needs  to  be  examined  relative  humidity  distribution  and  current  density  distribution  to  improve  the  durability. 
In  this  study,  the  PEFC  reaction  and  thermal  flow  analysis  model  including  the  effect  of  cooling  water  was  developed.  Furthermore,  the  effects 
of  the  shape  of  separator  channels  and  the  flow  pattern  on  current  density  distribution  and  relative  humidity  distribution  were  examined  by  this 
numerical  analysis.  As  a  result,  it  was  found  that  the  uniform  current  density  distribution  did  not  directly  relate  to  the  uniform  humidity  distribution, 
because  the  humidity  distribution  was  complexly  affected  by  the  generated  water,  the  water  transfer  between  the  anode  and  the  cathode  and  the 
gas  temperature.  Moreover,  it  was  confirmed  that  the  optimal  flow  pattern  of  gas  and  cooling  water  could  make  the  relative  humidity  higher  and 
more  uniform.  This  calculation  model  can  help  us  to  design  the  optimal  separator  shape  and  to  determine  the  optimal  operating  conditions  which 
can  narrow  the  low  humidity  area  and  improve  the  cell  durability. 

©  2006  Elsevier  B.V.  All  rights  reserved. 
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1.  Introduction 

Recently,  environmental  pollution  and  destruction  of  ecosys¬ 
tem  has  been  getting  worse  because  of  mass  consumption  of 
fossil  fuels  such  as  petroleum.  The  exhaustion  of  these  energy 
resources  becomes  a  serious  problem.  In  order  to  contribute 
to  solve  these  problems,  fuel  cells  are  expected  to  be  practical 
use  because  it  emits  less  environmental  pollutant  and  converts 
more  efficiently  from  chemical  energy  to  electrical  energy  than 
other  energy  resources.  Especially,  polymer  electrolyte  fuel  cell 
(PEFC)  is  expected  as  driving  power  of  vehicles  and  stationary 
power  supply,  because  it  can  work  at  low  temperature  and  has 
high  power  density.  The  performance  of  PEFC  has  improved 
rapidly  by  developing  the  new  component  materials  and  opti¬ 
mizing  the  system.  In  order  to  spread  PEFC  for  various  uses,  it 
is  necessary  to  improve  the  durability  and  the  cell  output,  and 
to  reduce  the  cost.  The  power  generated  by  PEFC  is  affected  by 
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the  structure,  the  material  and  the  operating  conditions  through 
the  process  of  generation.  The  phenomena  of  mass  transfer,  heat 
transfer,  catalysis,  electrochemical  reactions  and  fluid  dynamics 
are  shown  only  in  an  internal  cell,  and  it  is  greatly  important 
to  understand  the  correlation  among  such  complex  phenomena 
in  detail  to  improve  and  optimize  the  PEFC  component  and  the 
system.  However,  these  phenomena  are  caused  regardless  of  the 
size  of  the  area  from  an  interface  of  catalyst  layers  to  a  stack, 
and  these  phenomena  affect  each  other  intricately.  Therefore,  it 
is  very  difficult  to  measure  local  conditions  accurately  by  exper¬ 
iments,  and  very  few  researchers  examine  that. 

In  recent  years,  a  numerical  analysis  method  has  been  used 
to  examine  that.  Nguyen  and  White  [1],  and  Yi  and  Nguyen 
[2]  developed  the  heat  and  water  transport  models  (2D)  which 
accounted  for  various  operating  conditions  and  hydration  of  the 
membrane.  On  the  other  hand,  it  is  thought  that  the  analysis 
with  the  computational  fluid  dynamics  (CFD)  is  important  to 
calculate  the  transport  phenomena  in  detail,  and  such  kinds 
of  studies  have  been  increasing  recently.  Um  et  al.  [3]  and 
Wang  et  al.  [4]  have  developed  the  two-dimensional  model  with 
CFD  which  included  two-phase  flow.  Dutta  et  al.  [5]  made  the 
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three-dimensional  computational  model  based  on  a  commer¬ 
cial  software  package,  Fluent.  Berning  et  al.  [6]  presented  the 
non-isothermal  and  three-dimensional  models,  and  calculated 
the  distribution  of  current  density  and  concentration  in  straight 
channels.  Um  and  Wang  [7]  compared  the  performance  of  the 
flow  in  straight  channels  with  that  in  interdigitated  channels  by 
the  three-dimensional  analysis.  Berning  and  Djilali  [8]  exam¬ 
ined  the  effect  of  porosity  and  thickness  of  gas  diffusion  layers 
in  the  straight  channels  by  the  three-dimensional  model.  These 
PEFC  numerical  analysis  models  contributed  to  the  optimization 
of  component  design  and  the  operating  conditions. 

PEFC  needs  to  be  improved  its  durability  in  order  to  commer¬ 
cialize  this  system.  There  are  some  reports  concerning  the  dura¬ 
bility  under  various  operating  conditions,  and  it  was  reported  that 
degradation  of  MEA  caused  the  cell  voltage  to  reduce  and  the 
electrolyte  membrane  to  brake.  Though  the  degradation  mech¬ 
anism  has  not  been  investigated  in  detail  yet,  it  is  thought  that 
the  attack  by  peroxide  radical  is  responsible  for  this  degradation 
phenomenon,  and  it  is  reported  that  this  degradation  is  acceler¬ 
ated  under  a  low  humidity  condition  [9-12].  On  the  basis  of  this 
experimental  result,  it  is  effective  to  design  the  optimal  operat¬ 
ing  conditions  and  shapes  which  can  narrow  the  low  humidity 
area  in  order  to  improve  the  PEFC  durability.  In  PEFC,  water 
and  heat  are  generated  by  the  electrochemical  reactions,  and  the 
relative  humidity,  which  is  a  function  of  the  local  temperature 
and  the  vapor  concentration,  is  not  uniform.  And  in  the  case  of 
a  stack,  the  effect  of  the  cooling  water  for  heat  recovery  has  to 
be  considered.  There  are  some  experimental  studies  concerning 
the  humidity  distribution  in  a  cell,  and  important  results  were 
reported  in  ref.  [13].  On  the  other  hand,  in  order  to  examine  the 
humidity  distribution  under  various  conditions  and  with  various 
shapes,  the  numerical  analysis  is  effective  from  the  viewpoints 
of  cost  and  time. 

In  our  past  researches  [14],  the  effects  of  changing  the  oper¬ 
ation  temperature,  the  humidify  temperature  and  the  hydrogen 
and  oxygen  concentration  in  the  supplied  gas  on  the  i-V  charac¬ 
teristic  of  a  small  PEFC  were  examined  experimentally.  For  the 
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Fig.  1.  Numerical  analysis  model  including  cooling  water. 


experiment,  we  developed  two  models:  one  was  the  PEFC  reac¬ 
tion  model  that  could  show  these  influences  on  PEFC  reaction 
characteristics;  and  the  other  was  the  PEFC  reaction  and  flow 
analysis  model  that  was  combined  with  the  thermal  flow  analy¬ 
sis.  With  this  PEFC  reaction  and  flow  analysis  model,  five  kinds 
of  separators  were  evaluated  from  the  viewpoints  as  follows:  gas 
flow  condition,  uniformity  of  current  density  and  temperature, 
reduction  of  pressure  drop  and  ejection  of  water.  In  ref.  [15], 
the  simplified  two-dimensional  PEFC  analysis  model  including 
flow  and  heat  transfer  of  cooling  water  was  made.  The  influence 
of  changing  the  thickness  of  the  membrane  and  the  GDL  on  the 
cell  performance  was  calculated.  However,  these  previous  mod¬ 
els  of  ours  did  not  include  the  effect  of  the  gas  flow  through  the 
gas  diffusion  layers.  In  refs.  [16,17],  the  mass  transfer  and  the 
flow  in  the  gas  diffusion  layers  were  calculated,  and  the  approx¬ 
imate  model  for  the  GDL  mass  transfer  based  on  the  theoretical 
model  was  developed.  Next,  with  this  model,  the  PEFC  reaction 
and  thermal  flow  analysis  model  which  enabled  us  to  calculate 
an  actual-sized  cell  was  made.  The  numerical  analysis  made  it 
possible  to  examine  how  the  separator  depth  and  the  GDL  effec¬ 
tive  porosity  and  the  GDL  permeability  and  the  flow  rate  of  the 
cathode  gas  effected  on  the  output  performance  and  the  current 
density  distribution.  In  this  study,  a  calculation  of  temperature 
of  cooling  water  was  combined  with  our  past  model,  and  PEFC 
reaction  and  flow  analysis  model  including  the  effect  of  cooling 
water  was  developed.  Furthermore,  the  relative  humidity  dis¬ 
tribution  and  current  density  distribution  were  examined  under 
various  conditions  with  this  model. 


Table  1 

Equations  for  calculation  of  current  density  and  oxygen  mass  transfer  in  GDL 
Current  density 

V  =  E-  (RT/at2F)  ln[!C«f/(o2C^]  -  Rohmi,  y  =  io2/Cg 

Oxygen  mass  transfer  in  GDL 
Upstream:  Sh  =  fi  +  XRe0  5Sc0  8 
Downstream:  Sh  =  j3  +  X(Re  —  co)0-5 Sc0-5 

Conditions:  Sft  =  ((/4FX/GDL/(Dgf(C'h2annel-C®2))),  Re  =  IqdlpU/ fi,  Sc  =  n/pD^,  A  =  l/2y/GDL/7rffGDL,  Df  =  sDj 
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2.  Numerical  analysis  model  including  the  effect  of 
cooling  water 

2.7.  Single  cell  model  and  stack  model  with  cooling  water 

Fig.  1  shows  a  single  cell  analysis  model  used  in  this  study. 
As  shown  in  this  figure,  the  gas  flow  velocity,  the  concentra¬ 
tion  and  the  temperature  were  calculated  in  the  gas  channels  on 
the  anode  and  the  cathode.  And  it  was  assumed  that  the  tem¬ 
perature  distribution  of  MEA  and  GDL  were  the  same  as  each 
other  and  they  were  unified,  and  the  temperature  and  the  cur¬ 
rent  density  were  calculated  in  a  unified  part.  (This  unified  part 
is  described  as  a  solid  phase  from  now  on.)  In  this  model,  the 
temperature  distribution  of  cooling  water  flowing  behind  a  sep¬ 
arator  was  calculated  with  the  heat  transfer  among  the  anode 
gas,  the  cathode  gas,  the  solid  phase  and  the  cooling  water. 
Fig.  2  shows  a  stack  and  cooling  water  model.  The  cooling 
water  flows  among  each  cell.  The  heat  of  overvoltage  is  con¬ 
ducted  to  the  cooling  water  through  the  anode  gas,  the  cathode 
gas  and  the  MEA.  The  cooling  water  is  heated  at  an  outlet,  and 
this  warmed  water  is  utilized  for  cogeneration  in  the  case  of 
a  stationary  PEFC  system.  In  order  to  calculate  the  heat  bal¬ 


ance  of  cooling  water,  the  data  of  the  temperature  distribution 
of  adjoining  cells  is  needed.  However,  it  is  difficult  to  calcu¬ 
late  the  temperature  distribution  of  all  cells  and  the  cooling 
water  in  a  stack  composed  of  about  100-200  cells,  so  the  inlet 
condition  and  the  internal  condition  are  assumed  to  be  equal 
to  the  adjoining  cells  with  each  other  in  this  study.  In  Fig.  2, 
the  cell  A  and  the  cell  B  are  the  same  as  each  other,  and  the 
heat  balance  between  the  cooling  water  and  a  cell  were  virtually 
treated  as  periodic  condition.  It  is  thought  that  this  hypothesis 
is  effective  at  the  center  of  a  stack  which  consists  of  many  lay¬ 
ered  cells,  and  this  study  examines  the  center  of  a  stack.  The 
periodic  condition  at  the  both  ends  of  a  stack  is  inappropriate, 
and  so  the  model  improving  heat  balance  is  needed  in  the  next 
study. 

2.2.  Hypotheses  and  basic  equations  of  this  study  model 

The  governing  equations  in  this  simulation  were  derived  from 
the  following  assumptions: 

1 .  The  gas  flow  rate  at  the  inlet  in  each  channel  is  uniform. 
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Fig.  2.  Stack  and  cooling  water  model  in  a  numerical  analysis. 
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2.  The  volume  of  liquid  water  is  ignored  and  the  water  moves 
with  gas. 

3.  The  reduction  of  the  reacting  areas  caused  by  flooding  of 
electrodes  is  ignored,  and  it  is  also  ignored  that  liquid  water 
prevents  the  diffusion. 

4.  Fluid  is  incompressible  Newtonian  fluid  and  ideal  gas.  The 
flow  condition  is  laminar  flow,  and  the  gas  properties  are 
constant. 

5 .  The  heat  transfer  between  separators  and  gas  is  ignored.  But 
the  heat  transfer  among  gas  phase,  solid  phase  and  cooling 
water  is  included. 

6.  The  cell  voltage  is  uniform  and  constant. 

7.  Only  the  resistance  overvoltage  and  the  water  transfer  in 
membrane  include  the  influence  of  temperature. 

8 .  In  membrane,  ionic  conductivity,  electroosmosis  coefficient 
and  water  effective  diffusion  coefficient  that  depend  on 
membrane  humidity  are  determined  by  the  water  activity 
in  an  anode  side. 

9.  The  crossover  gas  through  membrane  is  disregarded. 

10.  The  permeability  of  GDL  is  constant  and  uniform. 


In  this  study,  the  one-dimensional  analysis  as  plug  flow  in 
each  channel  was  available  on  the  assumption  that  the  distribu¬ 
tion  of  the  gas  flow  rate  at  the  inlet  was  uniform.  Though  the 
separator  shape  was  a  two-dimensional  structure  to  the  direc¬ 
tion  of  membrane  surfaces,  the  quasi-two-dimensional  analysis 
model  could  be  applied  on  the  definition  that  the  direction  from 
the  inlet  to  the  outlet  was  a  positive  v  direction  in  each  serpentine 
channel.  As  a  result,  the  equations  were  simplified  and  it  took 
less  time  to  calculate.  However,  in  the  case  of  the  calculations  of 
the  temperature  distribution  on  the  solid  phase,  it  was  calculated 
by  the  two-dimensional  analysis  model.  In  other  studies  [18,19], 
it  was  reported  that  supplied  gas  flowed  through  GDL  which  was 
porous  media,  and  that  the  distribution  of  gas  flow  rate  in  chan¬ 
nels  was  not  uniform.  In  our  past  studies  [16,17],  we  developed 
the  analysis  model  including  the  effect  of  gas  flow  through  GDL 
on  cell  performance.  In  order  to  calculate  the  gas  flow  rate  that 
flowed  to  the  next  channel  through  GDL,  pressure  distribution 
was  calculated  by  the  two-dimensional  analysis.  Moreover,  the 
following  two  terms  in  gas  channels  were  ignored  to  simplify 
the  calculations;  the  heat  conduction  term  in  the  energy  balance 
equations;  and  the  diffusion  term  in  the  mass  balance  equations. 
In  the  calculations,  as  the  relationship  between  relative  humidity 
and  the  decreasing  rate  of  cell  voltage  was  unknown,  the  inter¬ 
nal  condition  of  a  cell  was  calculated  without  degradation  of 
MEA. 

The  equation  of  continuity  is  shown  by  the  following 
equation, 


dvk 

dx 


where  v  is  the  velocity  of  mixed  gas,  v  the  distance  along  a  gas 
flow  channel,  7?rea  all  reaction  rate,  Qb  all  gas  flow  rate  through 
GDL  per  unit  volume  to  the  next  channel  and  the  superscript  k 
is  the  anode  side  or  the  cathode  side.  7?rea  is  calculated  by  the 


following  equation. 


where  l^g  is  the  depth  of  gas  channels,  p  the  density  of  mixed 
gas,  Mj  the  molecular  weight  of  the  chemical  species  j  and  rj  is 
the  reaction  or  condensation  rate  per  unit  area  of  the  chemical 
species  j.  Qb  is  calculated  by  the  following  equation, 

q\  =  EeU  (3) 

n 


where  Q\^  is  the  gas  flow  rate  through  GDL  to  the  n  direction. 
The  equation  of  motion  is  shown  by  the  following  equation, 


-v/  +  Pkvk(R/d  +  Ql) 


where  p  is  the  pressure,  /i  the  gas  viscosity,  wq  the  width  of  gas 
channel  and  the  operator  D/D t  is  the  substantial  time  derivative 
that  is  shown  by  the  following  equation. 

Dnk  dvk  k  dvk 

- =  —  +  uk —  (5) 

D  t  dt  dx 

The  viscous  term  in  Eq.  (4)  is  derived  from  the  Hele-Show 
model.  This  term  includes  the  effect  of  the  viscous  drag  between 
two  pairs  of  facing  walls  in  a  channel. 

The  equation  of  the  chemical  species  j  is  shown  by  the  fol¬ 
lowing  equation, 


D  t 


-  A  +  C'jiR^  +  of) 

ld,g 

+  ECA«)2L,in)-ECAL,out) 

n  n 


where  Cj  is  the  concentration  of  the  chemical  species  j ,  Cj(n) 
the  concentration  of  the  chemical  species  j  at  the  next  chan¬ 
nel  to  the  n  direction,  Qb(n, in)  the  gas  flow  rate  through  GDL 
from  the  n  direction  adjoining  channels  to  this  point  and 
Qb(n, out)  is  the  gas  flow  rate  through  GDL  from  this  point  to 
the  n  direction  adjoining  channels.  The  equations  of  the  fol¬ 
lowing  eight  chemical  species  are  derived:  C^2,  C^2,  Cjj  0^, 

ch2o(1)’  Co,-  cn2>  CH20(vr  ch2o(1)  that  are  hydrogen,  oxy- 
gen,  nitrogen,  vapor  and  condensed  water  in  anode  and  cathode 

channels. 

The  equations  of  energy  are  shown  by  the  following  equa¬ 
tions, 


(Gas) 


Drk 

~dT 


+  ch  |  TkIRk  . 

pkCkZk  +  T  ^rea  +  Qb> 

+ErnkeLm)-ErkeLou,) 

n 


n 
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(Solid) 


scs  _  frSylj'S  +  +  *?4  +  *?!  +  ^6 


The  reaction  and  condensation  rates  of  each  ingredient  are 
shown  by  the  following  equations, 


(Cooling  water) 


D  rw 
Dt 


In  the  energy  equation  of  gas,  Cp  is  the  specific  heat,  T  the 
temperature  and  q\  and  q2  are  the  heat  fluxes  from  the  solid 
phase  and  the  cooling  water,  respectively.  Tn  is  the  gas  tem¬ 
perature  in  the  next  channel  to  the  n  direction.  In  the  equation 
of  solid  phase,  k  is  the  heat  conductivity,  /s  the  thickness  of 
solid  phase,  <73  the  heat  value  per  unit  area  as  a  result  of  elec¬ 
trochemical  reactions,  q\  and  <75  the  heat  fluxes  from  gas  and 
cooling  water,  respectively,  q§  the  latent  heat  flux  of  water 
condensation  and  the  superscript  s  is  the  solid  phase.  In  the 
equation  of  cooling  water,  /w  is  the  depth  of  channel  in  which 
cooling  water  flows  and  q ™  and  q ™  are  the  heat  fluxes  from 
both  sides  of  gas  channels  and  from  solid  phase,  respectively. 
These  heat  fluxes  and  heating  values  are  shown  by  the  following 
equations, 


q\  =  hk(Ts  -  Tk) 

q\  =  U^(TW  -  Tk ) 

q*3  =  (Eah  ~  V)i 

q\  =  h\Ta  -  Ts)  +  hc(Tc  -  Ts ) 

qs5  =  2C/|(TW  -  Ts) 

^6  =  — AHh2o(^h20(1)  +  rH20(l)) 
q™  =  [/f  (T2  -  Tw )  +  U^(TC  -  Tw ) 
q%  =  2U\(T%  -  rw) 


(10) 


where  h  is  the  heat  transfer  coefficient  of  anode  or  cathode 
gas,  Eah  the  heating  value  converted  into  theoretical  voltage, 
V  the  operating  voltage,  i  the  local  current  density,  A//h2o 
the  evaporation  enthalpy  of  water,  Uj  the  overall  heat  trans¬ 
fer  coefficient  between  the  anode  gas  or  the  cathode  gas  and 
the  cooling  water  and  Uj  is  the  overall  heat  transfer  coeffi¬ 
cient  between  the  cooling  water  and  the  solid  phase.  These 
overall  heat  transfer  coefficients  are  shown  by  the  following 
equations, 


T  (l//zk)  +  (/seP/£seP)  +  (l//!w) 

|  (11) 

U  S  —  _ 

T  ((/SeP+4g)/^SeP)  +  (l//!w) 


where  /sep  is  the  thickness  of  the  separator  between  the  cooling 
water  and  the  gas  phase  and  ksep  is  the  heat  conductivity  of  the 
separators. 
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h2o(v)  -  _<a  +  2o^2F  +  ^’g^c  (  Ch2o(v) 
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where  F  is  the  Faraday’s  constant,  R  the  gas  constant,  FH2o,sat 
the  saturated  vapor  pressure,  a  the  water  transfer  coefficient  and 
bc  is  the  condensation  rate  constant. 

The  gas  flow  rate  through  GDL  is  calculated  by  the  following 
Darcy’s  model, 


Q 


k  _ 
b(n ) 


k  Jk 
lGDL 

P  ify/w 


k 

L 


(13) 


where  kp  is  the  permeability  of  GDL,  /qdl  the  thickness  of  GDL 
and  wl  is  the  width  of  the  land  area  between  channels. 

The  current  density  distribution  was  calculated  by  the  equa¬ 
tion  obtained  in  our  past  study  [16].  The  equations  to  calculate 
the  current  density  are  shown  in  Table  1 .  The  approximate  equa¬ 
tion  of  the  oxygen  mass  transfer  rate  to  the  electrode  through 
GDL  was  obtained  by  the  numerical  analysis  and  the  theoretical 
model.  The  equation  was  the  function  of  the  Reynolds  number 
and  the  Schmitt  number  and  was  in  proportion  to  the  square  root 
of  the  Reynolds  number.  In  the  equations  shown  in  Table  1,  it 
was  found  that  the  /3  was  0.624  and  co  was  1.3  in  the  case  that 
the  width  of  a  channel  and  the  land  were  both  1  mm  and  the 
thickness  of  GDL  was  300  p,m.  The  current  density  was  calcu¬ 
lated  by  the  oxygen  mass  transfer  model  and  the  overvoltage 
equations  which  were  substituted  for  the  following  three  local 
factors:  the  concentration  of  oxygen  and  hydrogen  and  vapor, 
the  temperature  of  membrane  and  anode  gas  and  cathode  gas, 
and  the  gas  flow  velocity  in  GDL.  The  resistance  of  the  pro¬ 
ton  transfer  through  electrolyte  membrane  was  calculated  by 
the  Nguyen’s  equation  [1].  In  the  calculation  of  current  density 
with  equations  shown  in  Table  1 ,  the  same  values  as  reference 
[16]  were  used. 

It  is  thought  that  water  moves  with  electroosmosis  and  back- 
diffusion  in  the  electrolyte  membrane.  When  one  proton  moves 
from  the  anode  side  to  the  cathode  side,  the  water  movement 
coefficient  a  shows  the  net  number  of  water  molecules  moving 
along  with  protons.  This  is  from  the  method  by  Nguyen  and 
White  [1]. 
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The  local  concentration,  the  temperature,  the  flow  veloc¬ 
ity,  the  current  density  were  calculated  with  Eqs.  (1)— (13)  and 
Table  1 .  The  partial  differential  equations  are  solved  by  the  finite 
differential  method.  The  boundary  conditions  of  the  flow  veloc¬ 
ity,  the  temperature  and  the  concentration  are  set  as  follows: 

(1)  The  inlet  boundary  of  gas:  these  variables  are  constant. 

(2)  The  outlet  boundary  of  gas:  the  gradients  of  these  variables 
are  constant. 

The  current  density  and  the  water  transfer  coefficient  were 
calculated  all  over  the  electrode  areas.  Those  variables  were 
calculated  until  they  became  stationary  state.  The  relative  errors 
of  the  balance  equations  of  mass,  species  and  energy  became 
less  than  1%  in  all  calculations.  In  this  study,  the  gas  flow  rate 
was  automatically  set  so  that  the  utilization  was  constant  to  the 
prescribed  average  current  density,  and  the  flow  rate  of  cooling 
water  was  automatically  set  so  that  the  temperature  at  the  outlet 
was  constant. 

3.  Results  and  discussions 

The  calculation  conditions  are  shown  in  Table  2.  Assum¬ 
ing  that  the  influence  of  the  changes  of  the  gas  composition  on 
the  physical  property  is  very  little,  the  physical  properties  were 
treated  to  be  constant.  The  shapes  of  the  gas  channels  are  shown 
in  Fig.  3.  The  electrode  area  is  225  cm2  (a  square,  15  cm  on  a 
side),  and  the  widths  of  the  channels  and  the  land  are  both  1  mm. 
A  is  an  ordinary  serpentine  separator,  B  is  a  distributed  serpen¬ 
tine  separator  and  both  separators  have  25  channels.  The  shape 
of  the  cooling  water  channels  is  the  same  as  that  of  anode  gas 
channels  or  cathode  gas  channels.  Fig.  4  shows  the  flow  pattern 
of  the  gas  and  the  cooling  water  and  Table  3  shows  the  combi¬ 
nation  of  separator  shape  and  flow  pattern.  In  an  actual  PEFC 
stack,  as  the  position  of  the  manifold  is  restricted,  it  is  difficult  to 
set  an  inlet  (or  an  outlet)  of  anode  and  cathode  at  the  same  posi- 


(A) 


Table  2 


Operation  condition  and  shape  of  cell 

Pressure  (MPa) 

0.1 

Inlet  gas  temperature  (°C) 

70 

Inlet  cooling  water  temperature  (°C) 

70 

Outlet  cooling  water  temperature  (°C) 

75 

Humidify  temperature  (°C) 

65 

Inlet  gas  composition 

Anode 

H2:N2  =  75:25 

Cathode 

02:N2  =  21:79 

H2  utilization  (%) 

70 

O2  utilization  (%) 

40 

Thickness  of  membrane  (fxm) 

30 

Size  of  catalyst  layer  (cm2) 

225 

GDL  thickness  (fxm) 

300 

Channel  width  (mm) 

1 

Channel  depth  (mm) 

1 

Land  width  (mm) 

1 

GDL  permeability  (m2) 

2.5  X  10“u 

Electromotive  force  (V) 

1.23 

Average  current  density  (Acm~2) 

0.45 

tion  with  each  other.  Accordingly,  these  positions  were  set  at  the 
symmetric  position  to  each  other  in  this  study.  The  simulations 
were  carried  out  under  16  conditions  about  the  flow  pattern  of 
the  gas  and  the  cooling  water,  and  it  took  4  h  per  one  calculation 
by  penitum4®  3.2  GHz  PC. 

Fig.  5  shows  the  cell  voltage  and  the  average  relative  humidity 
of  the  anode  and  the  cathode  on  each  condition.  In  this  figure,  it 
was  found  that  the  cell  voltage  was  almost  equal  to  that  in  each 
condition.  The  average  humidity  of  the  counter  gas  flow  (nos.  3, 
4,  7,  8,  11,  12,  15  and  16)  was  about  5%  higher  than  that  of  the 
parallel  gas  flow  (nos.  1, 2, 5, 6, 9, 10, 13  and  14)  in  both  cases  of 
the  anode  and  the  cathode.  Since  the  average  current  density  and 
the  amount  of  generated  water  were  constant,  it  was  thought  that 
this  difference  was  caused  by  the  current  density  distribution, 
the  water  balance  between  the  anode  and  the  cathode  and  the 
temperature  distribution. 


(B) 


Fig.  3.  Separator  shape:  (A)  is  the  ordinary  serpentine  separator  with  25  channels  and  (B)  is  the  distributed  serpentine  separator  with  25  channels. 
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Counter  flow 
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Counter  flow 
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Anode  gas 


Cathode  gas 
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Fig.  4.  (a-d)  The  flow  pattern  of  gas  and  cooling  water. 


Table  3 


Combination  of  separator  shape  and  flow  pattern 


Number 

Anode 

separator 

Cathode 

separator 

Flow 

pattern 

Cooling  water 
pattern 

1 

A 

A 

Parallel 

Anode 

2 

A 

A 

Parallel 

Cathode 

3 

A 

A 

Counter 

Anode 

4 

A 

A 

Counter 

Cathode 

5 

A 

B 

Parallel 

Anode 

6 

A 

B 

Parallel 

Cathode 

7 

A 

B 

Counter 

Anode 

8 

A 

B 

Counter 

Cathode 

9 

B 

A 

Parallel 

Anode 

10 

B 

A 

Parallel 

Cathode 

11 

B 

A 

Counter 

Anode 

12 

B 

A 

Counter 

Cathode 

13 

B 

B 

Parallel 

Anode 

14 

B 

B 

Parallel 

Cathode 

15 

B 

B 

Counter 

Anode 

16 

B 

B 

Counter 

Cathode 

Combination  No. 

Fig.  5.  Calculation  results  of  cell  voltage  and  average  relative  humidity. 
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Fig.  6.  Current  density  distribution  under  each  flow  condition. 


The  current  density  distribution  on  each  condition  is  shown 
in  Fig.  6.  In  this  figure,  the  distribution  of  the  parallel  gas  flow 
(nos.  1,  2,  5,  6,  9,  10,  13  and  14)  was  more  remarkable  than  that 
of  the  counter  gas  flow  (nos.  3,  4,  7,  8,  11,  12,  15  and  16).  It 
confirms  that  the  influence  of  the  flow  pattern  of  cooling  water 
on  the  current  density  distribution  was  not  large.  In  addition, 
it  was  found  that  the  difference  of  humidity  between  a  right 
side  of  the  distributed  serpentine  separator  and  a  left  side  of  it 
was  larger,  and  that  the  current  density  distribution  was  more 
uneven  when  the  distributed  serpentine  separator  was  used  in 
the  cathode  channels. 


The  anode  and  the  cathode  relative  humidity  distribution  are 
shown  in  Figs.  7  and  8,  respectively.  In  these  figures,  it  could 
be  confirmed  that  the  trend  of  the  humidity  distribution  in  both 
sides  were  almost  equal  to  each  other.  Therefore,  it  was  thought 
that  the  humidity  distribution  in  both  sides  were  related  to  each 
other  by  electroosmosis  and  back-diffusion,  and  that  one  did 
not  become  more  uneven  when  the  other  was  comparatively 
uniform.  In  addition,  the  humidity  distribution  of  the  anode  was 
more  remarkable  than  that  of  the  cathode  in  each  condition,  so 
it  was  thought  that  the  anode  humidity  distribution  affected  the 
cell  durability  more. 
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Fig.  7.  Anode  relative  humidity  distribution  under  each  flow  condition. 


Next,  the  influence  of  the  flow  pattern  of  cooling  water  on  the 
anode  relative  humidity  distribution  was  examined.  In  Fig.  7, 
when  both  gases  were  counter  flow  and  the  cooling  water  was 
synchronized  with  the  cathode  gas,  the  anode  humidity  distri¬ 
bution  was  uniform  comparatively.  The  reasons  are  thought  as 
follows: 

(1)  In  the  upstream  area  of  the  anode,  the  vapor  concentration 
increased  because  of  back-diffusion  from  the  downstream 
area  of  the  cathode  where  the  vapor  concentration  was  high. 


And  the  gas  temperature  rose  because  the  cooling  water  was 
synchronized  with  the  flow  of  the  cathode  gas. 

(2)  In  the  downstream  area  of  the  anode,  the  vapor  concentration 
decreased  because  of  electroosmosis  and  the  diffusion  to  the 
upstream  area  of  the  cathode  where  the  vapor  concentration 
was  low.  And  the  gas  temperature  fell  because  the  cooling 
water  was  synchronized  with  the  flow  of  the  cathode  gas. 

(3 )  The  anode  vapor  concentration  and  the  gas  temperature  were 
both  high  at  the  upstream  area  of  the  anode  (1)  and  these 
were  both  low  at  the  downstream  area  of  the  anode  (2).  Con- 


G.  Inoue  et  al.  /  Journal  of  Power  Sources  162  (2006)  81-93 


91 


Flow 

pattern 


Parallel 


Counter 


Cooling 

water 


Anode 


Cathode 


Anode 


Cathode 


Separator 

Anode:A 

Cathode:A 


Separator 

Anode:A 

Cathode:B 


7 


8 


fc^aassT— -  _ - |  -  lij 

ll^| 

Separator 

Anode:B 

Cathode:A 


Separator 

Anode:B 

Cathode:B 


13 


14 


4 


16 


Anode  gas 


70  80  90  100 

Cathode  relative  humidity  [%] 

— ►  Cathode  gas  — | 


Cooling  water 


Fig.  8.  Cathode  relative  humidity  distribution  under  each  flow  condition. 


sequently,  the  relative  humidity  was  uniformly  distributed 
all  over  the  electrode  area. 

It  is  also  possible  to  examine  the  cause  of  the  humidity  dis¬ 
tribution  by  investigating  the  water  balance  between  the  anode 
and  the  cathode.  Next,  the  influences  of  channel  shapes  were 
examined.  The  difference  of  humidity  between  a  right  side  of 
the  distributed  serpentine  separator  and  a  left  side  of  it  was  larger. 
The  unevenness  was  the  largest  in  particular  when  the  gas  flow 
pattern  was  parallel  and  the  flow  pattern  of  cooling  water  was 
synchronized  with  the  flow  of  the  anode  gas  (nos.  1,  5,  9  and 
13).  Whichever  the  combination  was,  the  relative  humidity  was 


lower  at  the  inlet  of  the  cathode,  and  that  was  higher  at  the  inlet 
of  the  anode  under  this  operating  condition.  Therefore,  the  bun¬ 
dle  of  the  gas  channels  was  restricted  to  the  left  side  (or  the  right 
side)  of  the  electrode  area  because  of  the  distributed  serpentine 
separator,  and  such  uneven  distribution  was  formed. 

In  order  to  investigate  the  current  density  and  humidity  dis¬ 
tribution  quantitatively,  the  standard  deviations  of  the  current 
density  and  humidity  were  calculated.  Fig.  9  shows  them  under 
each  condition.  In  this  figure,  it  was  confirmed  that  the  charac¬ 
ters  of  the  anode  and  the  cathode  humidity  under  each  condition 
were  the  same  as  each  other,  and  that  these  characters  were  not 
the  same  as  those  of  the  current  density.  It  was  thought  that 
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Fig.  9.  Standard  deviation  of  current  density  and  relative  humidity  under  each 
condition. 


the  uniform  current  density  distribution  did  not  directly  relate 
to  the  uniform  humidity  distribution,  and  that  it  resulted  from 
that  current  density  distribution  was  affected  by  the  oxygen  con¬ 
centration  and  the  ionic  conductivity  which  was  the  function  of 
humidity. 

In  Fig.  9,  the  standard  deviation  of  humidity  in  no.  4  is  the 
lowest.  Fig.  10  shows  the  component  ratio  of  the  anode  relative 
humidity  distribution  on  each  condition.  In  this  figure,  the  low 
humidity  area  (less  than  70%)  in  no.  4  was  smaller  than  that 
in  other  conditions.  Although  the  minimum  permissible  relative 
humidity  to  prevent  degradation  of  membrane  cannot  be  decided 
because  the  degradation  rate  has  not  been  studied,  it  is  expected 
that  the  combination  of  no.  4  can  make  the  relative  humidity 
more  uniform  and  higher,  and  that  this  combination  improves 
the  cell  durability.  In  addition,  as  the  relative  humidity  of  this 
combination  was  less  than  90%,  it  was  expected  that  the  effect  of 
the  inhibition  of  the  gas  flow  by  liquid  water  was  lower  than  that 
of  other  cases.  From  the  reasons  mentioned  above,  the  relative 
humidity  in  no.  4  is  the  highest  and  the  most  uniform. 
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Fig.  10.  Component  ratio  of  anode  relative  humidity  distribution. 


4.  Conclusion 

The  PEFC  reaction  and  flow  analysis  model  including  the 
effect  of  heat  management  by  cooling  water  applied  in  a  PEFC 
stack  was  developed.  The  optimal  separator  shape  and  the  opti¬ 
mal  flow  pattern  of  the  gas  and  the  cooling  water  that  make  the 
relative  humidity  higher  and  more  uniform  were  examined  with 
this  model  under  each  condition.  The  following  results  were 
obtained  by  these  examinations: 

1 .  When  the  gas  flow  was  counter,  the  average  relative  humidity 
was  larger  than  that  of  other  cases. 

2.  As  the  humidity  distribution  in  the  anode  and  the  cathode 
sides  were  related  to  each  other  by  electroosmosis  and  back- 
diffusion,  the  trend  of  the  humidity  distribution  in  both  sides 
were  almost  equal  to  each  other.  And  the  distribution  of  the 
anode  was  more  remarkable  than  that  of  the  cathode  in  each 
condition. 

3 .  The  uniform  current  density  distribution  do  not  directly  relate 
to  the  uniform  humidity  distribution.  It  resulted  from  that  the 
current  density  distribution  was  affected  by  the  oxygen  con¬ 
centration  and  the  ionic  conductivity  which  was  the  function 
of  humidity,  and  that  the  relative  humidity  was  affected  by 
the  generated  water  and  the  gas  temperature. 

4.  In  the  case  of  the  distributed  serpentine  separator,  the  differ¬ 
ence  of  humidity  between  the  right  side  and  the  left  side  was 
larger. 

Furthermore,  the  relative  humidity  is  the  highest  and  the 
most  uniform  in  the  following  cases;  that  the  gas  flow  pattern 
is  counter,  that  the  cooling  water  is  synchronized  with  the  flow 
of  the  cathode  gas,  and  that  the  ordinary  serpentine  separator 
is  used  in  the  anode  and  the  cathode  side.  Using  the  numerical 
analysis  model  developed  in  this  study,  it  is  possible  to  examine 
the  relative  the  humidity  distribution  that  affects  the  cell  dura¬ 
bility  and  to  evaluate  the  optimal  design  to  control  humidity 
distribution  under  various  operating  conditions. 
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